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A  Theory  For  The  Sraall  Uns^nmnetrlc  Deformations  of 

Cylindrical  Shells 

1,    Introduction 

In  reference  1  approximate  theories  to  determine  the 
rotationally  symmetric  deformations  of  "thick"  cylindrical  shells 
were  developed  in  a  systematic  manner  from  the  exact  three  dimen- 
sional theory  of  linear  elasticity  for  an  isotropic  and 
homogeneous  material  (hereafter  referred  to  as  the  exact  theory). 
These  approximate  theories  were  obtained  by  a  method  that  is  in 

a  sense,  a  generalisation  of  the  boundary  layer  expansion 
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technique  used  by  Priedrichs  and  later  Priedrichs  and  Dressier 

in  a  study  of  the  bonding  of  plates."   However  the  scaling  and 

expansion  techniques  differ  somewhat  from  those  in  reference 

and  3« 

In  this  paper,  the  method  presented  in  reference  1  is 

extended  to  include  unsymmetric  deformations  of  cylindrical  shells. 

Sirailar  scaling  and  expansion  procedures  are  employed  and  the 

present  results  specialize  to  those  of  reference  1  for  rotationally 

symmetric  deformations.   The  basic  element  in  the  approximation 

is  a  "thin"  shell  theory.   The  "thick"  shell  theories  are  then 

obtained  as  ''corrections"  to  the  thin  shell  theory,   u'e  find 

that  the  so  called  Donnell  theory  of  cylindrical  shells  ,  which 

appears  naturally  as  a  consequence  of  oior  expansion  procedures, 

is  the  appropriate  thin  shell  theory,  A  summary  of  the  thick 

shell  theory  is  given  in  Section  6. 

"In  this  connection  see  also  reference  Ij., 


^ 

2-4   Formulation, 

In  terms  of  the  cylindrical  coordinate  system, (r,  G,  x),  a 
cylindrical  shell  of  Liniform  thickness  .2h,  radius  R  and  length 
L  is  defined  as  an  elastic  body  bounded  by  the  coaxial  surfaces 
r  =  R  1  h  and  the  planes  x  =  0  and  x  =  L.   The  inner  and  outer 
surfaces  of  the  cylinder,  r  =  R  -  h  and  r  =  R  +  h,  respectively, 
are  subjected  only  to  normal  pressures.   The  edges  of  the  cylinder 
formed  by  the  annular  regions  x  =  0,  L,  R  -  h  <  r  <  R  +  h,  are 
submitted  to  normal  and  shear  forces.   The  applied  forces  are 
arbitrary  provided  they  form  a  system  in  equilibrium* 

Our  pxirpose  is  to  obtain,  in  a  systematic  manner  from  the 
exact  theory,  tvjo  dimensional  or  "shell''  theories  which  approx- 
iriiate  the  deformations  of  cylindrical  bodies  with  "small''  values 
of  h/R.   It  is  convenient  to  Introduce  the  dimensionless 
variables: 


X  _  X      A  =  ~  =    ^        ^  _  r-R  ^  r-'R 
R&^    /Rh  '      e^   (h/R)^/^  '      Rz 


where,  a,  is  a  positive  integer  and  e  is  the  dimensionless 
parameter, 

(2)  e=   1)1/2^  , 

The  inner  and  outer  surfaces  are  thus  given  by  ^  =  ]|^  1, 

respectively,  while  the  ends  of  the  cylinder  correspond  to 

C  =  0,  L'  where,  L'  =  — -  =  ,  A  more  general  scaling 

Re^    \/Rh 
procedure  than  (1)  and  (2)  may  be  used.   However,  this  leads  to 

"unreasonable"  shell  theories  vjhich  are  discussed  in  reference  1 
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for  the  special  case  of  rotationally  symmetric  deformations. 
Attention  in  this  paper  is  therefore  restricted  to  the  scaling  in 
(1)  and  (2),   The  integer,  a,  provides  some  flexibility  in  the  choice 
of  the  parameter,  e.  For  simplicity  of  presentation  we  set  a  =  1 
and  hence," 

&  =  (h/R)^/^  .    ' 

Dlmensionless  stress  and  displacement  components  are  defined 
by  dlvidinfj  the  physical  stresses  by  Young's  modulus  and  the  physical 
displacements  by  R.   Using  these  dlmensionless  quantities  and  the 
variables  (1)  and  (2),  the  equations  of  the  exact  theory  neglecting 
body  forces  are,  with  an  obvious  notation: 
equilibrium  equations, 

-1, 


■■^-   J  =  0  , 


(3) 


^zx,C-*-^^x,5^S^^^^^zx-*-^"^xM^  =°' 


+£-^<=^  ^]  -  0  ; 


stress  displacement  relations  (Hooke's  Law), 


^In  prin'clple,  no  difficulty  is  enco-untered  if,  a,  is  an  arbitrary 
integer.   However  the  calculations  become  more  cumbersome ^  see 
reference  1. 
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v,i    +  ew  =  ed  +  e^^)[<^  -  v{cr    +  <r  )]    , 


.   +  (1  +  c^^)v,^  =  e(l  +  e^J:)2(l  +  v)^^ 


v,^  -  g(^)[v  -  £"-'-v;,^]   =  £^2(1  +  v)cr^  ^ 
Compatibility  Equations,  • 

A^^  +  g(^)  /a,^  +2g^(^)(cr^  -   o^   +  2£-^a5Q^^ 
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where   v  is  Poisson's  ratio  and. 
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In  (3»5)  we  have  employed  the  conventional  notation  that  an 
independent  variable  appearing  as  a  subscript  following  a  coimna 
denotes  differentiation  ^^rith  respect  to  that  variable.  Thus, 
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A  complete  formulation  of  the  elastic  problem  for  the  cylinder 
(Formulation  A)  consists  of  (3)  and  (Li.)  and  appropriate  boundary 
conditions  on  the  siirfaces  and  ends  of  the  cylinder,  A  second 
and  equivalent  formulation  (Formulation  B)  consists  of  (3)  and  (5) 
with  appropriate  boundary  conditions.   These  conditions  may  be 
obtained  by  specifying  the  applied  forces  in  the  following  manner: 


(6a) 


f<r^yiA>  t  lj> 


£)  =  <3:^(5,4,  t  lie)  =  0  , 


^zo 


cA^AAi^)   =  Vr^^^Ai^),  cr(?*4>  -  i;e)   =  Pt^^*^*^)* 
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^(0,4, ^je)  =  ^(cl.,^;£),  o;^(0,4,^;e)  =^(f),^;e)  , 


<^^(0,4,^;£)  =  t((l.,^ie)  ; 


and  conditions  similar  to  (6b)  on  ^  =  L,   Here  p  (5,4>e-)  and 
p_(5,(})^£)  are,  respectively,  the  applied  normal  forces  on  the 
outer  and  inner  surfaces  of  the  cylinder,  ^|,^;e),  ^{i(,K,»^)   and 


'"  c(4j  -  l»e)  =  t(4i,  -  l^e)  =  0  for  continuity, 


t '? 


■i.ir>.^z.S,?>.- 


.vJ-£;If;i^Cfc;> 


'?^ -;i -iSiiifK-- 


:  a 


J--    -.. 


(3J.# 


V" 


^O^i 


l^':.-^^r.~■^ 


^:-^-^^-^^:  v-Ji?4,.o]:.;;> 


.    /  •■■•«•<•• 


v.'/-  ^(a^^vf):^    <^'^'^^^^^  ^^'*  "^^^  s 


=Tr: 


?A- 


r-  .?■  .V' 


"^((Ij^je)   Sive   the  forces  applied  to  the   ends   of  the   cylinder,     I7e 
further   impose   the   unessential  but   convenient  restriction  that. 


(7)  Po  ^  Pi  • 


Formulation  A  is  employed  in  the  following  sections  to  discuss 
the  deformations  in  the  "interior"  of  the  cylinder.   To  examine 
the  deformations  in  the  "boundary  laj'-er""  we  use  Formulation  B 
in  sections  l^   and  5« 

3»    The^  Interior  Problems.. 

V/e  assume  that  each  component  of  stress  and  displacement, 
indicated  by  the  generic  symbol,  O'i^t^ jKfS.) ,   can  be  represented, 
asymptotically,  as  a  pov;er  series  in  e: 

oo 

(8)  (T{E,A*Ki^)^IZo'''i^>^»K)^'' 

n=0  .   ■   ■   , 

where  a^{^,^,V   =  0  if  n  <  0,   The  functions  cr^{^,^,K)   are  called 
the  interior  stress  coefficients  or  the  interior  displacement 
coefficients  of  order  n  which  ever  the  case  may  be.  VJe  further 
assume  that  the  prescribed  forces  in  (6)  can  be  expanded  in  povrer 
series  in  ei 


The  term  "boundary  layer"  refers  to  a  narroxir  volume  of  the 
cylinder  adjacent  to  and  including'  the  ends  where  stresses  and 
deflections  may  change  rapidly  in  the  ^  direction.   The  remaining 
part  of  the  cylinder  is  called  the  "interior". 


.'i^i-; 


.TT. 


.to  XL! 


0   ,U«.n-!^^ii..p  I  .. 


rr>  y-  ".-  o        c,-  -j!^.'  T 


00 


^?««5tii.:^^ 


.:■  >  n  li  q  T  (;;.f. 


-.r  1, 


/  ■ 


(9. a) 


P^(C,4i£) 


CO  f>'^^ 


S  P]:(^,4;£)    n=o  Ip^C^,!)  I 


n 


(9.b) 


5^(4, ^*e) 


A 


'^^{^,k;^) 


> 


00 


v^ 


t(4,?;e) 


n=0 


-< 


nrn 


t^(c|i,^|£) 


Expansions  of  the  form  (8)  are  substituted  into  the  equil- 
ibrium equations  (3)   and  the  stress  displacement  relations  ([{.)• 
Coefficients  of  the  same  pox-iers  of  &  are  equated  yielding  a  system 
of  differential  equations  satisfied  by  the  interior  stress  and 
displacement  coefficients  of  all  orders.   The  coefficients  of  e^ 
are : 


.n 


(10a) 
(10b) 
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(lie)       v%   +  ^v^:2  +  u^  =  2(1  +  v)[<-^  +  K<^1^]      , 


(lid)        w?^=^^-^.  V (0-^-2  -.02"^)   , 

(lie)        u?^  +  w?:^  =  2(1  +  v)tf^;2   , 

(llf )        v?^  +  ?v5-2  .  v^-2  +  w^f  =  2(1  H-  v)[<T^-2  f  ^p^^-^   . 


Here, 


i+2(  j+l)^ 

i+2j+l=n         ' 

Appropriate  boiindary  conditions  on  ^  =  -  1  for  the  interior  stress 
and  displacement  coefficients  are  obtained  by  substituting  (9a) 
and  expansions  of  the  forn  (8)  into  (6a)  and  equating  coefficients 
of  like  powers  of  e.   This  yields,  for  the  coefficients  of  e^, 

(12a)  ^zx^^'^'  -  ^^  ='^zO^^'^'  1  1)  =  0  , 

(12b)         ^cr'li^s^A),  crli^f^^'DJ  ^^P^i^A),     p5(?,4)1 

Equations  to  determine  each  interior  coefficient  are  obtained 
from  (10-12).  Ue  first  solve  (11a, b,c)  f or  <r^,  ^j-^  andcr!^^, 
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If  n  =  0,    it   folloi-js   from  (11),   recalling  that  (^  =  0  i£  n  <   0, 
that. 


o  o  o     _     o     _  ,  o 


Therefore, 

(111)  u°(^,4,^)   =  v°(5,4,^)   =  0,w°(5,4,^)   =  1J°(5,4)      , 

where  W^('l,<Ji)  is  an  arbitrary  fimction  to  be  subsequently 
deterirdned.'"   Similarly  vje  find  frora  (10b, c)  and  (12a)  that, 

(i5a)  <rly,^^»^>^)  =crlQ^iAsO  =  0    . 

Using  this  result  in  conjunction  with  (12b)  and  (7),  Sq.  (10a) 
yields. 


"Here  we  have  assurr.ed  that  u°(^,4,'^)  is  periodic  in  (j>  and  have 
neglected  rigid  body  deformations. 
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(l5b) 


<^2  =  pS  =  P?  =  °'     ^^  ^  <  2- 


Returning  to  (lid)  we  see  that. 
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Employing  this  last  result,  (II4.)  and  (l5a),  (lle,f)  yields 
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uV.  +  i^V-"^  =  v,^  +  w,^  =  0, 
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which  implies  that. 
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(17) 


^      if  n  =  1,2  . 


v^(5,4,^)  =  -¥?^-^^  +  v^(^,4). 


Here  U  (ii,^)  and  V^(^,(|))  are  arbitrary  finictions  which  are 
determined  later. 

Substitution  of  (1^)  and  (17)  into  (13)  gives  the  first  two 
stress-displacement  relations  for  o^,   c^  and  O^q   in  terms  of 
U-'-,V-'-  and  17°,  and  U^,V^  and  W-'-j 
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VJe  now  obtain  two  systems  of  differential  equations  to 
determine  U-'-,V-'-  and  W°,  and  U^,V^  and  ¥■*•,   Substituting  (l5a)  and 
(18)  into  (10b, c)  i-jith  n  =  1  or  2,  integrating  vjitli  respect  to  ^ 
and  using  the  surface  conditions  (12a)  there  results. 
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The  coefficients  defined  in  (19)  are;  proportional  to  the  conven- 
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Equations  (20)  supply  four  of  the  required  differential  equations. 

Equations  (21)  are  the  stress  displacement  relations  for  the  first 

two  non-vanishing  interior  transverse  shear  stress  coefficients. 

The  remaining  tvjo  differential  equations  and  stress  displacement 

relations  for  the  first  two  non-vanishing  0~    are  obtained  by  sub- 

z 

stituting  (15),  (18)  and  (21)  into  (10a)  with  n  =  2  or  3. 
Integrating  the  resulting  equation  once  v/ith  respect  to  ^  and 
using  the  surface  conditions  (12b)  we  find  that. 
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+  s2~^^  +  (Pq  +  P^)  ,    if  n  =  2,3   . 

It  is  noi-J  assumed,  with  no  loss  of  generality,  that  both  p  (4,4) 
and  Pj(?,4)  ci-O  ^o"t  simultaneously  Identically  vanish.  This  implies 

that  the  normal  forces  applied  to  the  inner  and  outer  surfaces  are 

2   h 
of  order  of  magnitude  e  ~  r» 

To  summarize,  we  have  obtained  by  a  systematic  expansion 
procedure  relations  expressing  each  of  the  first  tvio  non-vanishing 
interior  stress  coefficients  in  terms  of  the  first  two  non- 
vanishing  interior  displacement  coefficients.   In  addition  approp- 
riate differential  equations  to  determine  tliese  displacement 
coefficients  have  been  obtained.   These  results  may  be  written  in 
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a  more  familiar  form.   Substituting  (19)  into  (20)  and  (2a)  and 
setting  n  =  0  gives  the  differential  equations  of  the  "zeroth 
order  interior  problem"  as: 


(2i^) 


lA  J      l-V     Trl  4.     1+V     ttI  4.    1  tO  _      n. 


A^w°  +  3(vi^  +  w°  +  vu^,)  =  3(1  -  v2)(p2  -  p2)     . 

while  from  (lO),  (19),  (21)  and  (23)  the  stress  displacement 
relations  of  the  zeroth  order  interior  problem  are. 
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^         2(l-v"^) 


+  (w?, .  +  vw;„,)(l  -  ^2) 


'4ct 


'^•■3 


+  2(v}^  +  1-J°  +  vU^^)^]  +  (p2  +  p2)   . 


J.O     c 


-  £L  ;i-:.;-if--:-i 


>-r  ■ 


;>^  V    i^,:.^'??  ^   !.^,K?^?^ 


■:■    {i,i)^^ 


r  1 .'  *^i'   J- 


-^'<w^>**.„-'r-^~' 


.,  sO-r^ 


I  c.  \   /^    -    ..-  '     '<   /  .>. 


/    %   I' — '-^ 


k'  ^l-ll 


ih  -  i 


V  +    ^^?Vi)   + 


'  1  - 


Ik. 


Here  /\     is  the   two-dimensional  biharraonio   operator 


and  the  coefficients  in  (2^)  are  given  in  (19).   The  differential 

equations  and  stress  displacement  relations  of  the  ''first  order 

interior  problem"  are  similarly  obtained  from  (18)  -  (23)»  These 

expressions  will  be  referred  to  as  Eq,  (27),   They  are  identical 

with  (2[j.)  -  (26)  provided  that  a,  1,  is  added  to  the  superscript 

of  each  term.   If  the  deformations  are  rotationally  symmetric  then 

all  the  previous  results  reduce  to  those  of  reference  1, 

The  differential  equations  (2I4.)  and  the  stress  displacement 

relations  (2^)  are,  in  our  notation,  identical  with  those  of  the 

thin  shell  theory  of  Donne  11,''  Equations  (26),  which  are  not  given 

in  the  Donnell  theory,  provide  "thick  shell  corrections"  to  this 

thin  shell  theory.   Thus  the  Donnell  equations  appear  in  a  natural 

way  as  part  of  the  first  approximation  of  the  exact  theory," 

The  differential  equations  and  stress  displacement  relations 

for   interior    problems   of   order  tvjo   or   greater  are   obtained,    in  a 

similar  fashion,  from  (10  -  13 )•  Because  of  the  complexity  of  the 

analysis  they  are  not  explicitly  shovm  here.  We  merely  remark  that, 

p   2   2    3      3  2 

<^,  (TQi  o~-^Qf   "i-^  ^nd  v  are  cubics  in  ^,  while  w  is  quadratic  in  ^. 

To  complete  the  formulation  of  the  interior  problems,  boundary 

conditions  appropriate  to  the  differential  equations  (20)  and  (22) 

and  (27)  are  required.   These  boundary  conditions  are  systeraatically 


^This  does  not  necessarily  imply  that  other  thin  shell  theories,^"* 
could  not  serve  as  first  approximations  to  the  exact  equations  if 
different  expansion  procedures  or  parameters  are  employed. 


oj.ict   bxjtIT      ,     iw 


7     II 
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derived  in  sections  I|.  and  ^  from  the  exact  theory.  V'e  observe  that 
all  the  results  for  the  interior  problems  are  obtained  vjithout 
reference  to  the  edge  boundary  conditions,  (6b )» 

i|»   Formulation  of  the  Boundar:/'  Layer  Problem. 

The  specific  dependence  of  the  interior  stress  coefficients 
<r"',  cr^  and  (7%^  on  ^  is  given  by  the  Interior  analysis."   However 
the  boundary  conditions  on  the  ends  of  the  shell,  (6b),  imply  that 
these  coefficients  are  arbitrary  functions  of  ^»   Therefore,  the 
expansions  (8)  cannot  satisfy  these  conditions  and  if  they 
represent  the  solution  they  do  so  only  in  a  region  away  from  the 
ends,  i.e.  (8)  are  not  uniformly  valid.   To  obtain  expansions  that 
are  valid  up  to  and  including  the  end  c^  =  0,  we  introduce  a 

"boimdary  layer  coordinate", n ,  by  "stretching"  the  variable  ^   so 

I 
that  the  resulting  expansions  may  represent  the  solution  'unlf  ormlyj  ' 


(28)  ^  =  |  =  E  • 


Thus  by  making  e  sufficiently  small,  every  neighborhood  of  the  end 
in  the  ^  variable  corresponds  to  an  tirbitrarily  large  one  in  the 
Yi    variable.   The  (j)  variable  is  not  stretched.   This  implies  that 
the  rapid  variations  that  occur  near  the  boundary  talce  place  only 
in  the  t,   direction,  i.e.,  the  direction  normal  to  the  boundary. 

"""For  example,  see  (25)  and  (26). 

'For  descriptions  of  these  boundary  layer  methods  see  References 
1  -  Ij.'.   In  the  following  we  confine  our  attention  to  the  end 
4=0,   An  identical  analysis  is  valid  for  the  end  ^  =  L'. 


li/  ^,  jj      fiAi'i  J  ^ 


*x^ 


■Cf     ,   .^   5"9.t 


•      *  . 


•SJ-s     yJ 


xi 


F^r 


r  r ......    If  r-^. 


O    Bf'f 


*  \  C 
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We   introduce    (28)    into  formulation  B   of  Section  2  and  define 
the  boundary  layer   stresses   indicated  by  the   generic   symbol, 
f(n,(|,^je),   as   the   physical   stresses   divided  by  Yoijng's  modulus. 
Each  of  the   boundary  layer   stresses    is   assumed  to  be   representable 
by  a  power  series   in  e: 


(29) 


CO 


fir^A,Ki&)-^IZ  f^{r^A,Vi 


where  f"^(n,({),^)  are  the  boundary  layer  stress  coefficients  and 


a  _ 


0  if  i  <  0,   Substituting  (28)  and  (29)  into  the  equilibrlxm 


and  compatibility  equations  (3)  and  (5)  and  equating  coefficients 


n. 


of  the  same  powers  of  e,  x^^e  obtain  for  the  coefficients  of  e  : 


-.n 


(f^  -  4)  + 


'tif.1 

-  ^zQ   =  0 


(30)    A  ^x,S  *  f",„  *  I='4  ^ 


? 


-xO        ^     ' 


^Q,^  -^  ^xQ,r^-*-2ll'fL^  +IZ:"4  =  0   ' 


^^' 


■2^f>  "■    ~-».^>-  * 
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>i 


v^ie  +  r?..+  zzi"i  =  0 


■X        '    '7>^ 


v^i^  -  v\^  -  rz'fl  -  2  z 


-"'  -  .i 


(f;  - 


ffi)  +  i^ 


zQ~ 


f;^  =  o. 


■»r^i 


v^fS  +  iz  tfS■^  r"]  +2IZ  (f,  -  fj)  +i^iz"^f,a=o, 


■IV^l 


(31) 


V^f^x  ^  ^7^  ^  2='^ 


zx 


0 


V^r^.  ^  EZ'^o  -  n"  r  ^,  -  IZ"r^.  ^  IZ^^r^,  =  0    , 


v'C^ii:'^.^^ 


ri.- 1+ 


ZO 


+  z:^^i2(fj  -  f^)  -  r"]  =  0 . 


2  t 

Here,   the   binomial  expansions   of  g(^)   and  g   (^)   are   employed,  2. 

and  >         are   defined  following  Eq's.    (11),   and 


V     i,    J    >   0   , 


-IV^i 


i+2j+3=n  ^ 


2  _  a^  ^  5^ 


V^^ 


pn  _  _1    r-p^  J.  -pn    ,    ^r.n-| 


T    -. 


'        ,,>,•>■ 


-rtrxTXc  J: 


tj:"l':!-[. /::^>i  X 


^^:^(x  -h  u^{-»:"^"*  '3^ 
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To  complete  the  formulation,  the  boundary  layer  stresses  must 
satisfy  the  boundary  conditions  at  the  ends  of  the  cylinder  and 
on  the  inner  and  outer  surfaces.   In  addition  these  stresses 
should  approach  or  ''match"  the  interior  stresses  for  "large'' 
values  of  V,,  i.e.,  "small"  values  of  e. 

From  (6),  (9)  and  (29)  I'/e  find  that  on  the  boundaries  of  the 
cylinder  the  f^[yig^,i:l,)    satisfy  the  following  conditions: 


(32) 


(f2(v^,4,i),r^(^,^,.i)]  =  lJ|pS-^p?-^-}.^ 


(33) 


Jf5J(o,4,^)  =^i^sK)  s   fS^(o,4,^)  =r(c^,^)  , 


i^^{oA,a  =  t^(ct.,^)   . 


"IT  ■'- 

Here  the  functions  Pq,.!^)  and  Pj.(4)  are  given  by'": 


Poj(^)  = 


i 


1  ^Po(0>4) 


^C" 


0 


?Ui)  = 


]_.  dPj(0j4) 

Ti 


a€^ 


0 


j  =  0,1,..., 

j  =  0,1,..., 

J  <  0  . 


'"See  (9a). 


s^.n 
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They  are  the  coefficients  obtained  by  expanding  p  and  Pj  in  a 
Taylor  series  in  ^  about  ^  =  0,  The  formulas,  (32),  follow  by 
substituting  (28)  in  this  series  and  then  using  (29)« 

The  ''matching  conditions",  or  the  asymptotic  form  of  the 
boundary  layer  stress  coefficients,  for  a  fixed  <j),  are  written  a: 


Oij.) 


lim  f"(n,cl.,^)  =  c^{yiA,K) 
r^ — >oo   ^ 


n  -  0,1, • • e , 


where. 


'-•'>„■■ 


(35a) 


o^(^,<l',^)  =IZs^-J(4,?)r^^   , 


j=0 


and 


(35b) 


,[     bY^sK)   =     ^ 


i  >  0  , 


i  <  0  , 


The  matching  conditions  are  derived  by  assuming  that  in  the 
neighborhood  of  ^  =  0  each  interior  stress  coefficient  can  be 
expanded  in  a  Taylor  series  in  $.   The  conditions  (3^4-)  and  (35) 
follow"  by  substituting  (28)  into  this  expansion  and  using  (8) 
and  (29). 


''See  Ref,  1  for  details  of  the  derivation. 


.n 


l^iKr' 


i^^^i) 


(■-. 


0 


h 


^  i^.^}^^ 


;  b\ 


.■?/"il903  saGT.-J'a    'ioi-iOoryx  rlope"  0  -  J  "lo   faoo-u 
: ,   itolBHBqxe  aiti-;^  cctxii   {8S)  gdx-iydiJ-scfc'a- 

•  V  S.  ^i  / 


•-xioirfsvJ.'tf^t'  «irf  lo  ei 


5»   Analysis  of  the  Boundary  Layer  problein. 

In  analyzing  the  boundarj'-  layer  problem,  Eq«s«  (30  -  35),  we 
frequently  use  certain  integral  relations  deduced  from  the 
solutions  of  Dirichlet  problems  for  the  harmonic  and  biharraonio 
equations  on  the  semi-infinite  strip,  D  5  1^1  <  1*  oo  >  ^7  >  0, 
These  relations  are  given  in  the  Appendix,         ..  o    ■-  ,-  ■ 

For  every  value  of  n  equations  (30  -  35)  separate  into  two 

distinct  systems  which  we  denote  as  Problem  P  and  Problem  T, 

Problem  P  involves   only  ffj,  f^,  f^  and  f^^,  while  Problem  T 

X    z    w      zx 

considers  the  remaining  two  stress  coefficients*   For  n  =  0  we  have 
for  Problem  P  the  first  two  of  (30)  and  the  first  four  of  (31)^ 


(36) 


^sK 


+  f 


zx. 


-   0 


^zx,^  *   ^x,ip 


=  0 


(37) 


t2^o 


-< 


V^f 2  =  0   , 


2^o 


"^^x"^  r?nn  =  0 ,      v'f^  +  n 


^K 


=    0 


v'f°  *  r;,.x  =  0 . 


zx 


'V^ 


The  appropriate  boundary  and  matching  conditions  follow  from 
(32  -  35): 


(38) 


C(^,4,  ii)  =  t^i^A,  ID  =  0  , 


©v.' 


8B 


•  ^  ^  =i  ,v 


^'•y' 
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(39) 


[s2(o,4)  +  s°(o,(|.)^],[o],[o]?  . 


In  the  boundary  value  problem  (36  -  39)  the  variable  (j)  is  a 
parameter.   This  problem  can  be  reduced  to  one  for  the  biharmonic 
equation  on  D  by  introducing  the  reduced  boimdary  layer  stress 
coefficients,  P  ('7,4*^)*  a-nd  an  associated  Airy  stress  fiinction, 
yi^{v),^fK.]  f    in  the  follox/ing  manners 


F°  =  f  ° 
zx         zx 


By  direct  substitution  vie  see  that  (I|.0)  is  a  solution  of  (39  -  i+S) 
provided  that  X  is  a  solution  for  a  fixed  value  of  4  of  the 
following  boundary  value  problem  on  D: 


v^^°  - 


0 


ihl) 


lim   [%?r,^(^,4,^),  X?^«.(7,c}.,^)]  =  0  , 

yC,^^ioA,a  =^i^,^)  -  s°(o,4)  -  s°(o,(t)^  , 


t  V.,-. 


■;■  i'^  'li 


»-J.S5s)  C 


♦      ■:'<.   -■  » 


•c.tlol 


^v,/t;^N^</X 


. }   cfjeri*    .:sa   &trr  ncid-jjrt: 


;.-,i.o3   B  3X 


;C[  no  neldotq   eidisv  Y'^-'"' 


0  =^  'X'^V 


=^  (X  t  *Kv>-iytx  =  (I  t  .s'^ty-yptc 


0  =  [(2*Kv0v--?X  .(.^.J-tS;>)..va 

'      tJ«r>  ■■■■  /      C^  I 


00  < — 0 


.  sOr*  \0,-,        .'V    1  >  O'-  _    .  •,  \     c. 


X 


-s^' 


(Ili) 


.'V      l-.'^""  ...  ^Y 
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Application  of  integral  relations  I  (see  Appendix)  to  (Lp.)  yields, 

1  1 

(I^2a,b)     S^(0,(1,J  =|J  ^°((i),^)d^  ,    S°(0,(|j)  -I  \  C^°(t,^)c3^ 

"■'■  -1 


1 


Equations  (I|.2a,b)  are  the  first  t\<ro   of  the  required  four  bOLindary 
conditions  for  the  differential  equations  of  the  zeroth  order 
interior  problem,  (2i{.),   Equation  ([|.2c)  implies  that  X°{^,K>)   is  a 
self -equilibrating  force.   For  convenience  we  hereafter  exclude 
all  self-equilibrating  edge  forces  and  v/ithout  loss  of  generality 
assume  that  'cri^tK,)  d   0  and  '5'^ (4,^)  =   0. 

Problem  T,  for  n  =  0,  is  obtained  from  the  remaining  equili- 
brium and  compatibility  equations,  boundary  conditions  and  matching 
conditions  in  (30  -  3b)  vjith  n  =  0:        ■  ■■  "j.  :" 


Considering  4  as  a  parameter  the  problem  given  by  (I|.3)  can  be 
converted  into  Neum.ann's  problem  for  the  harmonic  equation  on  the 
semi-infinite  strip.   This  is  done  by  again  Introducing  reduced 
boundary  layer  stress  coefficients  and  an  appropriate  stress 
function  ^1/°(V7,4*^)*  such  that 


abXc. 


>    ■■ 


)    :    ^aoi..!-/5l«'x   S.£ 


:'l' 


,-.)''*' a 


ebwXi: 


'I/SS-f--. 


■>v-t  ftn/^ 


•aiid  no'X'. 


uj-do 


^     ft     '■ '  ,_ 


*nl:rfo^ 


y-k.'-         TV-V     ,- 


r  : 


::^  •^  .    -1 


^     U,|>)^B'  «    (^^KOJ.'^'i  ^    0   -    (I  t    .K 


^. 


<*  ■.  '-<• 


:s(<Ko)^::^-  ^-  i.Ko;^^..i,iu]|  :^   ^^..i,^ 


OiTii^' 
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Equation  {h^k)    is  a  solution  of  (^3)  if  \l/°  is  a  solution  of  the 
follov;ing  boundary  value  problem  on  D: 

M   t-  :  ^ ■   , ,  ;  , 

^^ioA,K.)   =  t°((t,^)  -  s°^(o,(|)  -  s°^(o,(l,)^  , 

ik^f  ^l'?^(7,4,  1  1)  =  0  ,  .   .  .  ^ 

lim  [^°        ^%]    =   0   , 
1^ — >oo   ^*         \ 

Existence  of  a  solution  of  {1^.$)   requires  the  integral  of  the  normal 
derivative  of  ^     aroiuid  the  boundary  to  vanish.   This  yields, 

1 

-1 

which  is  the  third  bovmdary  condition  to  be  associated  with  (21^), 
This  completes  the  analysis  of  (30  -  35)  for  n  =  0, 

V/ith  n  =  1  we  obtain  from  (30  -  35),  for  Problem  P,  equations 
and  boundary  conditions  identical,  except  for  superscripts,  with  ' 
(36  -  38).   The  matching  conditions  are  from  (34),  (35),  (l5b), 
(18)  and  (21), 


'Instead  of  using  the  stress  function  vj/  we  could  have  employed  a 
complex  conjugate  function.   This  leads  to  a  Dirichlet  problem 
for  the  harmonic  equation. 


.  "*  ox^  -  i^x-  ''     .'■■    "  g!  .  ■■"  '^'^' 


-J.v-'  ,       i; 


isXiJo' 


=    Cp?*    .;s;*!. 


I 


r 


•0  =  n  101   {ilZ.  -  0£.)  lo  aia^Isns  erl:^  aectelqrfioc 

-^-  >"  ■     ■    -ol   ^iZi   -  0£)  not?-  .C  =  n  r(:Hw 

,0.-k:I)    t(5C}    t(-l^-^J  ^-c>'il   T-.iij  snoidxijnoo  3nhioct.3iu  eriT     •(&£  -   r-M 

{.  (  X ; )    5n.£ 


•  i 
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(ij.7a) 

[sj(0,<|i)  +3^(0, (|))^  +  (S^^^COj)  +  s2^.(o,(l.)^)^  ],[0]  , 

VJe  acain  introduce  reduced  boundary  layer  stress  coefficients, 
F  (»o,cj),^)  and  an  associated  stress  function  '^i'^t^/O    through  the 
relations: 


^x  =  ^x  -  r^ji   =  ^^^  ^  '^^'4^ 


(47b) 


pi  =  f J  -  lira  fl   =  vV^X^  +  2(1  +  v)rl/?.   .  "    ' 
^ — >oo  ^ 

Equations  (I4.7)  provide  a  solution  of  the  equations  and  boiindary 
and  matching  conditions  of  Problem  P  with  n  =  1  if  ^  is  a  solution, 
for  every  4*  of  the  following  boundary  value  problem  on  the 
semi-infinite  strip; 


'       .'         r         »■ 


■    I       i       r  / 

.     -v    '•      -""^j     ^-y     *-'!■     ••  -  L        ,- f  (• 


1    1'0}J  ^(.^(K 


d.':..J    X"  . 


•■T-    ts 


"•  0-' 
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v^ — >oo    (^  ^^ 

Employing  (Ii.5),  the  third  of  the  Integral  relations  I  (see 
Appendix)  applied  to  (Ij.8)  yields  the  fourth  and  final  boundary 
condition  for  the  zeroth  order  interior  problem,  (2l|.)s 

1  1 

(1^.9)   f[2s°^,,t,(0*^)  +  s°^,(0,4)]  =^^^{^,K)dK.  +  \^  Kt°,^i^yK)ciK     . 

The  boimdary  conditions,  {l^2a,h),    (14.6)  and  (49)  coincide  with  those 
usually  associated  v;ith  the  Donnell  equations.   Equations  (l+Sa) 
and  (ii.6)  correspond  to  the  prescription  of  the  resultant  axial  and 
twist  forces  respectively,  V/hile  (14.2b)  and  (Ij.9)  are  respectively 
the  boundary  conditions  for  the  resultant  bending  moment  and 
"shear  reaction," 

The  remaining  two  integral  relations  applied  to  (I|.8)  give  the 
following  boundary  conditions  for  the  first  order  interior 
differential  equations  (27), 


'■-     ^i      i;    '  V;    •• 


'AK  c'- 


A' 


■iaiv-lx*^:^^ 


"ioi'i&cfri. 


'x^^ 


[(^^'^^    '"^  '   <: 


(^4^) 


■■■^    r  ri.i  ^    ^  j^.  ^^    ^  i.i ;    ,  f  -f ,  : 


'jfirn/o 


oj   iSitocvJj^ 


■sd  o\ 
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'i^i'  ,  'i"-'  -■  !i- 


v(  Wax  .i 
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1  1 

(50)  1  _         "i 

-1(0,4.)  =\K  K^{^sK)(yK   +  3  \  ^vl/?^(0,(}),^)d^  . 

The  remaining  equations  and  conditions  contained  in  (30  -  35) 
for  n  =  1  yield  the  follox-iing  Problem  T:  > 


^1     _   ^o      \-j2Jl      _      r-O 


^Z«i,^  "^  ^XP,17    "^0,4'    "^  "^XO    "'  '(|r^'  ...  - 


(51)   f^o^^'^'  -  1)  =  0  ,  fxe^°'^'^^  =  ti(4,^)  ; 

r   '. 

v^ii::oo{^^^'''«}={^^^'e,4^0'^)  -^s^,^^(o,i))(i-^2) 


3  , 


[S^^(0,<|)  +  S^^(0,(|.)^  +  (s°^^^(o,(i,)  +  s°Q^^(o,<}))^)7]  j. 

The  boundary  value  problem  given  by  (5l)  can  be  reduced  to  one  for 
Poisson's  equation  on  the  semi-infinite  strip  by  defining  reduced 
boiondary  layer  stress  coefficients  and  a  stress  function,  3!  (*t»4»^^* 
such  that, 

(52) 


•      '.••-■.v'.' 


i-    C; 


vv. 


*  :i6* 


IS    \i-\i//    J   -      v.-  V. ; 


'<*. 


JOO*^-^"^ 


;  f  ^(i»'^j)-, 


v-'-u 
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Here   we   have   used   (ij.0)    and   (ij.1).      Equations   (52)   give   a   solution 
of    (5l)    if  "7"      is  a   solution  of  the    following  boundary  value 
problem  on  D: 

00 


(53) 


V^  3I  =  (1  -  V)  r  V^7§4.^d^?     , 


•rl.l    ^.^     .     .1    /n    I^    ^  •?! 


Z7^(o*4,^)   =  -t-'(cl),^)   +  S^q(0,(|.)    +  S;^^(0,(|)^     ,    , 

where  we  have  employed  (l.|.l)  and  the  previous  assumption  that 
'C°((j>^^)  =  0,  Application  of  the  integral  relation  II  of  the 
Appendix  to  (53)  gives  the  third  boundary  condition  for  the 
first  order  interior  problem  asi  't 

iM  S^^(0,4)  =1  ^  t^(|,?)d^   .      --,.•... 

This  completes  the  analysis  of  (30  -  35)  for  n  =  1, 

The  remaining  boundary  condition  for  the  first  order  interior 
problem  and  boundary  value  problems  for  determining  the  boundary 
layer  stress  coefficients  of  order  tv:o  are  obtained  from  an 
analysis  of  (30  -  35)  with  n  =  2,   This  analysis  is  lengthy  and 
somewhat  involved  and  is  therefore  not  presented*  V/e  merely  list 
the  result  of  applying  the  third  of  the  integral  relations  I.  This 
gives  the  fourth  and  final  boundary  condition  for  the  first  order 
interior  problem  as: 
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JL  J. 

(55)     f[2sJo  .  (0,(1.)  +3^^^(0,(1.)]  ={^^{^,K)dK   +  V  t^4(c|>,^)d^  , 

Higher  order  approximations  may  be  obtained  by  examining 
(30  -  35)  with  n  >  2»  ■•..■. 

6,   A  Summary  of  the  Shell  Theory. 

In  the  previous  sections  we  have  obtained  by  a  systematic 
expansion  procedure  in  the  "small"  parameter  e,  a  sequence  of 
boundary  value  problems  whose  solutions  approximate  the  three 
dimensional  state  of  stress  in  a  cylindrical  shell.   The  approxi- 
mation can  be  carried  to  any  desired  degree.   The  basic  element 
of  the  approximation  is  a  thin  shell  theory  which  forms  part  of 
the  zeroth  order  interior  problem.   The  proper  boundary  conditions 
for  this  theory  are  obtained  from  the  boundary  layer  analysis.   It 
is  implied  by  the  expansion  procedure  that  the  higher  order 
approximations  serve  as  "corrections''  to  the  basic  thin  shell 
approximation  yielding  "thick  shell"  theories  of  increasing 
degrees  of  accuracy,  V/e  therefore  define  the  thick  shell  theory 
of  order  N  through  the  relations:  \ 


(56a) 


+  -|F^(?/e,cj,,^),P^(?/e,(l),^),F^^(S/e,(l.,^)j]s^  , 


( 56b  )         ^ N 
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(56c)      <^/^  =n^z(^'^'^^^''  '■HlK'^^/^A^K)     . 

Here,  cr^'   are  the  stresses  associated  with  the  N''  order  thick 
shell  theory,  and  cr    and  P  are  respectively  the  interior  and  the 
reduced  boiondary  layer  stress  coefficients  of  i   order. 

In  siunmary  we  list  below  the  bo'undary  value  problems  frora 
which  the  stresses  of  the  zeroth  and  first  order  thick  shell 
theory  may  be  computed. 

Zeroth  order  interior  problem; 

Differential  equations:   Equations  {2[j.)j 

Stress  displacement  relations;   Equations  (2^),  (26)  and  (19)j 

Boundary  conditions,  on  .^  =  0  (siinilar  ones  apply  to 
5  =  L»):  (i|2a,b),  ([j.6)  and  (q.9). 

The  differential  equations  of  the  zeroth  order  interior 
problem  can  be  formulated  in  a  somevihat  different  but  equivalent 
manner  by  introducing  an  ^.iry  stress  function  G  [c,,^)    such  that. 


(^7)  S°  =  G°^^,    Sl  =   gJ^^,    S°^  -  -G°^^   . 


Employing  (57)  v;e  find  that  the  first  two  of  (2[j.)  are  identically 

1      1 
satisfied,  see  also  (20),   Eliminating  U  and  V  from  the  first 

three  of  (19a)  with  n  =  0  we  find,  vising  (57),  that 


(58a)  AV  -  W?.^  =  0   , 
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and  that  the  third  of  (21;.)  becomes". 


(58b)        J^d°  +   3(i-v^)g;,^  =  3(l-v^)(p^  -  p|) 


The  differential  equations  {21^)    can  therefore  be  replaced  by  those 
in  (58)»  Appropriate  boundary  conditions  for  (58a)  are  obtained 
from  (57),  (I^2a)  and  ([j.6). 

Zeroth  order  boundary  layer  problems; 
Boundary  value  problems:   Equations  (ill )  and  ( 1-1-5  )| 
Stress  relations:   Equations  (39),  (i|0),  ([[.3)  and  ihj.\.) » 
First  order  interior  problem; 

The  differential  equations  and  the  stress  displacement 
relations  are  identical  with  those  of  the  zeroth  order  interior 
problem  provided  a,  1,  is  added  to  the  superscript  of  each  term. 
The  boundary  conditions  on  <|  =  0  are  given  by  (50),  (5^4-)  and  {$$) » 
Similar  conditions  apply  to  ^  =  L'. 

First  order  boundary  layer  problems: 

Boundary  value  problems:   Equations  (I)-8)  and  (53)* 

Stress  relations:   Equations  (1-1-7),  (5l)  and  (52)» 


'By  Introducing  the  complex  fimction,  lf{t,,^)   =  W°(^,4) 

+  iNy3(l-v  )C'  (4>4)*  the  system  (59)  may  be  combined  into  the 
single  complex  equation, 

A^H°  -  In/Ju^h;.,  =  3(l-v2)(p2  -  p2). 
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In  applying  the  thick  shell  theory  to  specific  probleins  it 
is  important  to  observe  the  order  of  solution.   The  zeroth  order 
interior  problem  must  be  solved  first.   The  solution  thus  obtained 
yields  Inhomogeneous  terms  for  the  zeroth  order  boundary  layer 
problems.   The  solutions  of  these  problems  in  turn  provides 
inhomogeneous  terms  for  the  first  order  interior  problem  and  so 
forth. 


-  HJ 
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Apnendix 

Consider  the  following  bovmdary  value  problem  for  X(f^,^) 
defined  on  D: 

V^X  =  0  , 

TC,^^(0,^)  =  f(^),  X,   ^^0,^)   =   g(^),   X,^^(<?,  +  1)  =  0  , 

lim  iX,^^'J^,K),     )c,t^^iyi,a]   =  0, 
7c,>^^(V,i)  =  k^(7). 

It  is  easy  to  show  that  if  X*y>^*  X'*?/  and  X»x  ^^®  continuous 
functions  and  if  Xf  Xfn  ^^'^  -X»y  ^^''^  single  valued  functions  then, 

1  09 

1  oo 

^1  0 


1 

i 


C  g(^)d^  =  0   , 


Consider  the  following  boundary  value   problem  for  _2_   i'ntK) 
defined  on  the    semi-infinite    strip: 


^! 


3k 


t2 


Z*^(o.^)  =  fU)  >        Z*J^*  i  1)  =  0  , 

It  Is  easy  to  show  tloat  if  3!»y  ^^  ^  continuous  function  and  if 
"-"  is  a  single  valued  function  then, 

1 

II  C  f (^)d^  =  0  . 
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